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Abstract 

We present cosmological perturbation theory in neutrinos probe interacting dark-energy models, 
and calculate cosmic microwave background anisotropics and matter power spectrum. In these 
models, the evolution of the mass of neutrinos is determined by the quintessence scalar field, which 
is responsible for the cosmic acceleration today. We consider several types of scalar field potentials 
and put constraints on the coupling parameter between neutrinos and dark energy. Assuming the 
flatness of the universe, the constraint we can derive from the current observation is ^ m^j < O.STeT^ 
at the 95 % confidence level for the sum over three species of neutrinos. We also discuss on the 
stability issue of the our model and on the impact of the scattering term in Boltzmann equation 
from the mass- varying neutrinos. 
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I. INTRODUCTION 

After SNIa{l| and WMApQ] observations during last decade, the discovery of the acceler- 
ated expansion of the universe is a major challenge of particle physics and cosmology. There 
are currently three candidates for the Dark-Energy which derives this accelerated expansion: 

• a non-zero cosmo logical constant {sl, 



• a dynamical cosmological constant (Quintessence scalar field) [4], 

• modifications of Einstein Theory of Gravity Sj 

The scalar field model like quintessence is a simple model with time dependent w, which 
is generally larger than —1. Because the different w leads to a different expansion history 
of the universe, the geometrical measurements of cosmic expansion through observations of 
SNIa, CMB, and Baryon Acoustic Oscillations (BAO) can give us tight constraints on w. 
One of the interesting way to study the scalar field dark energy models is to investigate the 
coupling between the dark energy and the other matter fields. In fact, a number of models 
which realize the interaction between dark energy and dark matter, or even visible matters, 
have been p.po.ed =o f„ fl B. Q . Ob™.o,. of ..e effects of .he.e ,n.e.c.o. 
will offer an unique opportunity to detect a cosmological scalar field 0, 

In this paper, after reviewing shortly the main idea of the three possible candidates of dark 
energy and their cosmological phenomena in section II, we discuss the interacting dark energy 
model, paying particular attention to the interacting mechanism between dark-energy with a 
hot dark-matter (neutrinos) in section III. In this so-called Mass- Varying Neutinos (MaVaNs) 



model 121], we calculate explicitly Cosmic Microwave Background(CMB) radiation and Large 
Scale Structure(LSS) within cosmological perturbation theory. The evolution of the mass of 
neutrinos is determined by the quintessence scalar filed, which is responsible for the cosmic 
acceleration today. Recently, perturbation equations for this class of models are nicely 



presented by Brookfield et al. 



131], (see also ij]) which are necessary to compute CMB and 



LSS spectra. A main difference here from their works is that we correctly take into account 
the scattering term in the geodesic equation of neutrinos, which was omitted there (see, 
however, [3]). We will show that this leads significant differences in the resultant spectra and 
hence the different observational constraints. In section IV, we discuss three different types of 
quintessence potential, namely, an inverse power law potential, a supergravity potential, and 



an exponential type potential. By computing CMB and LSS spectra with these quintessential 
potentials and comparing them to the latest observations, the constraints on the present mass 
of neutrinos and coupling parameters are derived. In conclusion we discuss two important 
points of this work on the impact of the scattering term of the Boltzmann equation and 
on the stability issue in the interacting dark-energy model. In appendix A, the explicit 
calculation for the consistency check of our calculations in section III is shown. Since we 
were asked to show explicit derivation of geodesic equation after our first draft was released, 
we show them in appendix B. 



II. THREE POSSIBLE SOLUTIONS FOR ACCELERATING UNIVERSE: 

Recent observations with Supernova la type (SNIa) and CMB radiation have provided 
strong evidence that we live now in an accelerating and almost fiat universe. In general, 
one believes that the dominance of a dark-energy component with negative pressure in the 
present era is responsible for the universe's accelerated expansion. However there are three 
possible solutions to explain the accelerating universe. The Einstein Equation in General 
Relativity is given by the following form: 

C^z. = = SttGT^^ + At/^^, (1) 

Here, G/j^i, term contains the information of geometrical structure, the energy-momentum 
tensor T^j, keeps the information of matter distributions, and the last term is so called the 
cosmological constant which contain the information of non-zero vacuum energy. After solve 
the Einstein equation, one can drive a simple relation: 

R AttG , , A 

5 = -— (p + 3p) + -. (2) 

In order to get the accelerating expansion, either cosmological constant A {loa — P/p — —1) 
becomes positive or a new concept of dark-energy with the negative pressure (cj^ < —1/3) 
needs to be introduced. Another solution can be given by the modification of geometrical 
structure which can provide a repulsive source of gravitational force. In this case, the attrac- 
tive gravitational force term is dominant in early stage of universe, however at later time near 
the present era, repulsive term become important and drives universe to be expanded with 
an acceleration. Also we can consider extra-energy density contributions from bulk space in 
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Brane- World scenario models, which can modify the Friedmann equation as oc p + p . In 
summary, we have three different solutions for the accelerating expansion of our universe as 
mentioned in the introduction. Probing for the origin of accelerating universe is the most 
important and challenged problem in high energy physics and cosmology now. The detail 



16| 



explanation and many references are in a useful review on dark energy 

In this paper, we concentrate on the second solution using the quintessence field. In 
present epoch, the potential term becomes important than kinetic term, which can easily 
explain the negative pressure with y9 —1. However there are many different versions 
of quintessence field: K-essence 17, 18|, phantom 19], quint omj2o|. ....etc., and to justify 



the origin of dark-energy from experimental observations is really a difficult job. Present 
updated value of the equation of states(EoS) are u = —1.02 ± 0.12 without any supernova 



data 



2l|. 



III. INTERACTING DARK-ENERGY WITH NEUTRINOS: 



As explained in previous section, it is really difficult to probe the origin of dark-energy 
when the dark-energy doesn't interact with other matters at all. Here we investi gate the 
cosmological implication of an idea of the dark-energy interacting with neutrinos 12|, |22 |. 
For simplicity, we consider the case that dark-energy and neutrinos are coupled such that 
the mass of the neutrinos is a function of the scalar field which drives the late time accel- 
erated expansion of the universe. In previous works by Fardon et al. 22 1 and R. Peccei 1^ . 
kinetic energy term was ignored and potential term was treated as a dynamical cosmology 
constant, which can be applicable for the dynamics near present epoch. However the kinetic 
contributions become important to descreibe cosmological perturbations in early stage of 
universe, which is fully considered in our analysis. 
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A. Cosmological perturbations: background Equations 

Equations for quintessence scalar field are given by 

(p + 2710 + a — = , (3) 

dcp 

Kff(0) = V{<P) + Vi{<P) , (4) 

VM) = a-' I -^,V¥+^^^)f{q) , (5) 

p <t> 

m,^(0) = rfiie *^pi (as an example), (6) 
where is the potential of quintessence scalar field, Vi{4>) is additional potential due to 



the coupling to neutrino particles [22, l23|, and m,y{(f)) is the mass of neutrino coupled to the 
scalar field. H is ^, where the dot represents the derivative with respect to the conformal 
time r. 

Energy densities of mass varying neutrino (MVN) and quintessence scalar field are de- 
scribed as 



-4 

Pu = a 



"^'^ Vq' + a^rnlfoiq) , (7) 



(27r)3 

,2 



3P, = a-' I -^_J^=/o(g) , (8) 
+ ^(0) , (9) 



= - V{<P) . (10) 

From equations ([7]) and ([8]), the equation of motion for the background energy density of 
neutrinos is given by 

r) In TD 

p. + snip. + Pu) = ^T-^0(p. - 3P.) . (11) 



B. Perturbation equations 

1. perturbations in the metric 

We work in the synchronous gauge and line element is 

ds'^ = a^{T) [-dr^ + {6ij + hij)dx'dx^] , (12) 



5 



In this metric the Chirstoffel symbols which have non-zero values are 



^ 00 



jk 



a 

a 

d ^ d 1 • 
-Oij + -hij + -a 
a a 



a 1 , 

a 
1 



(13) 
(14) 
(15) 
(16) 



where dot denotes conformal time derivative. For CMB anistropies we mainly consider the 
scalar type perturbations. We introduce two scalar fields, h{k,r) and r]{k,T), in k-space 



and write the scalar mode of hij as a Fourier integral 

hij{x,T) = 



3Q| 



kikjh{k, t) + {kikj - -5ij)6ry(fc, r) 



(17) 



where k = kk with A;*L = 1. 



2. perturbations in quintessence 

The equation of quintessence scalar field is given by 

□0- Kfr(0) = . (18) 

Let us write the scalar field as a sum of background value and perturbations around it, 
(f){x, t) = (I){t) + 50(x, r). The perturbation equation is then described as 

^50 + - ^V2((50) + + + 5 f ^1 = , (19) 



2a? 



where 



dVi _4 [ d^q de{q,(j)) 
e(g,</)) = v/g' + a2m2(0) , (21) 

d(p e{q,(p) dcj) ' ^ ' 

To describe 5 ^^j, we shall write the distribution function of neutrinos with background 
distribution and perturbation around it as 

/(a;\r,g,n,) = /o(r,g)(l + v[/(x\r,g,n,)) . (23) 



Then we can write 



dVi\ _4 f d?q d'^e _4 f d?q de 



where 



- 



\d(j) J 

For numerical purpose it is useful to rewrite the equations (12U1) and as 

(70 

Note that perturbation fluid variables in mass varying neutrinos are given by 



J (27r)3 e 7 (27r)3e2 

The energy momentum tensor of quintessence is given by 



and its perturbation is 



If .dV 



2 

This gives perturbations of quintessence in fulid variables as 

5p^ = -6T^ = \<j)6(f) + , 

d(j) 

SP^ = -5T°/3 = ^<P64> - 



d 



S} = T] - b]Tll'i = . 



C. Boltzmann Equation for Mass Varying Neutrino 



We have to consider Boltzmann equation to solve the evolution of VMN. A distribution 
function is written in terms of time (r), positions (x*) and their conjugate momentum [Pi). 
The conjugete momentum is defined as spatial parts of the 4-momentum with lower indices, 
i.e., Pi = mUi, where Ui = dxi/ {—ds^y^"^ . We also introduce locally orthonormal coordinate 
= (t, r*), and we write the energy and the momentum in this coordinate as {E,p^), w^here 
E = a/p^ + ml. The relations of these variables in synchronous gauge are given by [30], 



Po = -aE , 

Pi = a{6ij + ^hij)p' . 

Next we define comoving energy and momentum (e, qi) as 



(36) 
(37) 



e = aE = \/q'^ + a^ml , 
Qi = api . 



(38) 
(39) 



Hereafter, we shall use {x\q,nj,T) as phase space variables, replacing f{x\Pj,T) by 
f{x'^,q,nj,T). Here we have splitted the comoving momentum qj into its magnitude and 
direction: qj = qrij, where n^rii = 1. The Boltzmann equation is 



Df df ^ dx^ df ^ dqdf ^ drii df 



Dt dr dr (9x* dr dq dr drii 



dl 
dr 



(40) 



c 



in terms of these variables. From the time component of geodesic equation 25|, 



1 d 
2d^ 



{P^Y = -T%P''P^-mg"' 



m , 



(41) 



and the relation P^ 



.-2 



\/ q^ ^ a'^ml, we have 



dq !>• i j 
— = --hijqn n/ 
dr 2 



, m dm dx^ 
q dx^ dr 



(42) 



Our analytic formulas in eqs. (14111421) are different from those of [13| and IJ], since they have 
omitted the contribution of the varying neutrino mass term. We shall show later this term 
also give an important contribution in the first order perturbation of the Boltzman equation. 
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We will write down each term up to 0{h): 



df dfo ^ dfo , 

— = — + /o h — * 

dr dr dr dr 

dx^df_ _ q^i^ J 

dr dx^ e dx'^ ' 

dqdf ( ^^.^ydm^dx" 1 ■ .A 9/o 
= \ —a TT——, T^haqn w x — — 



We note that ^ and g are 0{h). 
1. Background equations 

Prom the equations above, the zeroth-order Boltzmann equation is 

f = 0. (44) 



The Fermi- Dirac distribution 



can be a solution. Here Qs is the number of spin degrees of freedom, h-p and k-Q are the 
Planck and the Boltzmann constants. 



2. perturbation equations 

The first-order Boltzmann equation is 

_iV-fcW!^^^|^^ = 0. (46) 

e q^ d(p omq 

Following previous studies, we shall assume that the initial momentum dependence is axially 
symmetric so that ^ depends on q — qh only through q and k • h. With this assumption, 
we expand the perturbation of distribution function, ^, in a Legendre series. 



Then we obtain the hierarchy for MVN 



where 



= H- t.3J^ , (48) 

e b omq 

^1 = liA;($o-2^2) + «:, (49) 



*, = -lM2*.-3*3)-^^ft+5^j g^, (50) 



« = . (52) 

3 e d(f) omq 



Here we used the recursion relation 

{i + l)P,+i(/i) = (2£ + l)/iP,(/^) - m^iif^) . (53) 
We have to solve these equations with a g-grid for every wavenumber k. 



IV. QUINTESSENCE POTENTIALS 

To determine the evolution of scalar field which couples to neutrinos, we should specify 
the potential of the scalar field. A variety of quintessence effective potentials can be found 
in the literature. In the present paper we examine three type of quintessential potentials. 
First we analyze what is a frequently invoked form for the effective potential of the tracker 
field, i.e., an inverse power law such as originally analyzed by Ratra and Peebles 



V^(0) = M^+"0-° (Model I) , (54) 



where M and a are parameters. 



We will also consider a modified form of V{(f)) as proposed by 39] based on the condition 



that the quintessence fields be part of supergravity models. The potential now becomes 

V{(P) = M^+"0-"e'*'/''"pi (Model II) , (55) 

where the exponential correction becomes important near the present time as — mpi. The 
fact that this potential has a minimum for = A/a/3mpi changes the dynamics. It causes 
the present value of w to evolve to a cosmological constant much quicker than for the bare 
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scale factor 



FIG. 1: Examples of the evolution of energy density in quintessence and the background fields 
as indicated. Model parameters taken to plot this figure are a = 10, 10, 1 for model I, II, III, 
respectively. The other parameters for the dark energy are fixed so that the energy densities in 
three types of dark energy should be the same at present. 

power-law potential [4^. In these models the parameter M is fixed by the condition that 
~ 0.7 at present. 

We will also analyze another class of tracking potential, namely, the potential of expo- 
nential type [4h] : 

V{<\)) = M^e"""^ (Model III) , (56) 

This type of potential can lead to accelerating expansion provided that a < \f2. In figure 
([1]) , we present examples of evolution of energy densities with these three types of potentials 
with vanishing coupling strength to neutrinos. 

A. Time evolution of neutrino mass and energy density in scalar field 

For an illustration we also plot examples of evolution of energy densities for interacting 
case with inverse power law potential (Model I) in Fig. ([2]). In interacting dark energy cases, 
the evolution of the scalar field is determined both by its own potential and interacting term 
from neutrinos. When neutrinos are highly relativistic, the interaction term can be expressed 
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FIG. 2: Examples of the evolution of energy density in quintessence and the background fields 
in coupled cases with inverse power law potential (Model I). Model parameters taken to plot this 
figure are a = l,/3 = l,3as indicated. The other parameters for the dark energy are fixed so that 
the energy densities in three types of dark energy should be the same at present. 



as 



10 
7^ 



(^QjTTli/ ^ /^^^massloss ' 



(57) 



where Pu^^^^^icss denotes the energy density of neutrinos with no mass. The term roughly 
scales as oc a~^, and therefore, it dominates deep in the radiation dominated era. However, 
because the motion of the scalar field driven by this interaction term is almost suppressed 
by the friction term, — 2?i0. The scalar field satisfies the slow roll condition similar to the 
infiation models, — 27i0 ^ a?^^^{p^ — 3P,y). Thus, the energy density in scalar field and the 
mass of neutrinos is frozen there. These behaviors are clearly seen in Figs. ([2]) and ([3]). 



B. Constrains on the MaVaNs parameters 



As was shown in the previous sections, the coupling between cosmological neutrinos and 
dark energy quintessence could modify the CMB and matter power spectra significantly. 
It is therefore possible and also important to put constraints on coupling parameters from 



current observations. For this purpose, we use the WMAP3 



35| and 2dF [36| data sets. 



The fiux power spectrum of the Lyman-a forest can be used to measure the matter 
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1000 




FIG. 3: Examples of the time evolution of neutrino mass in power law potential models (Model I) 
with Q = 1 and /? = (black solid line), (3 = 1 (red dashed line), (3 = 2 (blue dash-dotted line), 
(3 = 2> (dash-dot-dotted line). The larger coupling parameter leads to the larger mass in the early 



universe. 
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FIG. 4: The CMB angular power spectra for Model I. The solid line is the best fit for the model 
((a,/3) = (2.97,0.170)), the other lines are models with different parameter value of a and (3 as 
indicated. The points are WMAP three year data. 
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FIG. 5: The CMB angular power spectra for Model III. The solid line is the best fit for the model 
((a,/3) = (0.78,0.28)), the other lines are models with different parameter value of a and /3 as 
indicated. The points are 2dF data. 



power spectrum at small scales around z <3 42, |43|]. It has been shown, however, that the 
resultant constraint on neutrino mass can vary significantly from ^m^^ < 0.2eV to 0.4eV 
depending on the specific Lyman-a analysis used |44]. The complication arises because the 
result suffers from the systematic uncertainty regarding to the model for the intergalactic 
physical effects, i.e., damping wings, ionizing radiation fluctuations, galactic winds, and so 
on 451]. Therefore, we conservatively omit the Lyman-a forest data from our analysis. 



Because there are many other cosmological parameters than the MaVaNu parameters, 
we follow the Markov Chain Monte Carlo(MCMC) global fit approach js^] to explore the 
likelihood space and marginalize over the nuisance parameters to obtain the constraint on 
parameter(s) we are interested in. Our parameter space consists of 



P = {Qbh'^, flch'^, H, T, As, ris, rrii, a, (3) 



(58) 



where uibh"^ and VL(.h? are the baryon and CDM densities in units of critical density, H is 
the bubble parameter, r is the optical depth of Compton scattering to the last scattering 
surface, Ag and rtg are the amplitude and spectral index of primordial density fluctuations, 
and {mi, a, (3) are the parameters of MaVaNs deflned in section III. We have put priors on 
MaVaNs parameters as a > 0, and /5 > for simplicity and saving the computational time. 
Our results are shown in Figs. ([6]) - (IHI)- In these flgures we do not observe the strong 
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FIG. 6: Contours of constant relative probabilities in two dimensional parameter planes for inverse 
power law models. Lines correspond to 68% and 95.4% confidence limits. 
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FIG. 7: Same as Fig.®, but for SUGRA type models. 

degeneracy between the introduced parameters. This is why one can put tight constraints on 
MaVaNs parameters from observations. For both models we consider, larger a leads larger 
w at present. Therefore large a is not allowed due to the same reason that larger w is not 
allowed from the current observations. 
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FIG. 8: Same as Fig. ([6]), but for exponential type models. 

On the other hand, larger (3 will generally lead larger m^, in the early universe. This means 
that the effect of neutrinos on the density fluctuation of matter becomes larger leading to 
the larger damping of the power at small scales. A complication arise because the mass 
of neutrinos at the transition from the ultra-relativistic regime to the non-relativistic one 
is not a monotonic function of (3 as shown in Fig. ([3]). Even so, the coupled neutrinos give 
larger decrement of small scale power, and therefore one can limit the coupling parameter 
from the large scale structure data. 

One may wonder why we can get such a tight constraint on (3, because it is naively 
expected that large /3 value should be allowed ii Vtyh"^ ~ 0. In fact, a goodness of fit 
is still satisfactory with large j3 value when VtyK^ ~ 0. However, the parameters which 
give us the best goodness of fit does not mean the most likely parameters in general. In 
our parametrization, the accepted total volume by MCMC in the parameter space where 
Vtyh"^ ~ and {3 >1 was small, meaning that the probability of such a parameter set is low. 

We find no observational signature which favors the coupling between MaVaNs and 
quintessence scalar field, and obtain the upper limit on the coupling parameter as 

(3 < 1.11, 1.36, 1.53 , (59) 
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and the present mass of neutrinos is also limited to 

^^^Cday < 0.0095, 0.0090, 0.0084 , (60) 

for models I, II and III, respectively. When we apply the relation between the total sum of 
the neutrino masses My and their contributions to the energy density of the universe: VL^K^ = 
My/ (93.14eV^), we obtain the constraint on the total neutrino mass: My < 0.87eV{95%C.L.) 
in the neutrino probe dark-energy model. The total neutrino mass contributions in the power 
spectrum is shown in Fig[9l where we can see the significant deviation from observation data 
in the case of large neutrino masses. 

TABLE I: Global analysis data within la deviation for different types of the quintessence potential. 



Quantites 


Model I 


Model II 


Model III 


WMAP-3 data (ACDM) 




2.21 ±0.07 


2.22 ±0.07 


2.21 ±0.07 


2.23 ±0.07 




11.10 ±0.62 


11.10 ±0.65 


11.10 ±0.63 


12.8 ±0.8 


Ho 


65.97 ±3.61 


65.37 ±3.41 


65.61 ±3.26 


72 ±8 




10.87 ±2.58 


10.89 ± 2.62 


11.07 ±2.44 




a 


< 2.63 


< 7.78 


< 0.92 






< 0.46 


< 0.47 


< 0.58 




Us 


0.95 ±0.02 


0.95 ±0.02 


0.95 ± 0.02 


0.958 ±0.016 




20.66 ± 1.31 


20.69 ± 1.32 


20.72 ± 1.24 






68.54 ±4.81 


67.90 ± 4.47 


68.22 ±4.17 


71.6 ±5.5 


Age/Gyrs 


13.95 ± 0.20 


13.97 ±0.19 


13.69 ±0.19 


13.73 ±0.16 


nMVN h^[w^] 


< 0.44 


< 0.48 


< 0.48 


< 1.97(95%C.L.) 


T 


0.08 ± 0.03 


0.08 ± 0.03 


0.09 ± 0.03 


0.089 ± 0.030 



V. SUMMARY AND CONCLUSION 

Before concluding the paper we should comment two important points of this paper: the 
impact of the scattering term of the Boltzmann Equation in Sec. Ill and on the stability 
issue in the present models. 
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FIG. 9: Examples of the total neutrino mass contributions in power spectrum with Mi, 
0.9ey(Left panel) and with Mi, = 0.3ey (Right panel). Here the variable A is equal to a. 



Recently, 
field et al. 



perturbation equations for the MaVaNs models were nicely presented by Brook- 



13 1, (see also [14 1) which are necessary to compute CMB and LSS spectra. A 



main difference here from their works is that we correctly take into account the scattering 
term in the geodesic equation of neutrinos, which was omitted there (see, however, 15|). 



Because the term is proportional to ^ and first order quantity in perturbation, our results 
and those of earlier works |l^ remain the same in the background evolutions. However, 
as will be shown in the appendix, neglecting this term violates the energy momentum con- 
servation law at linear level leading to the anomalously large ISW effect. Because the term 
becomes important when neutrinos become massive, the late time ISW is mainly affected 
through the interaction between dark energy and neutrinos. Consequently, the differences 
show up at large angular scales. In Fig. ffTOl) . the differences are shown with and without 
the scattering term. The early ISW can also be affected by this term to some extent in some 
massive neutrino models and the height of the first acoustic peak could be changed. How- 
ever, the position of the peaks stays almost unchanged because the background expansion 
histories are the same. 



As shown in 
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47( 1 . some class of models with mass varying neutrinos suffers from the 



adiabatic instability at the first order perturbation level. This is caused by an additional force 
on neutrinos mediated by the quintessence scalar field and occurs when its effective mass is 
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FIG. 10: Differences between tfie CMB power spectra witfi and without the scattering term in the 
geodesic equation of neutrinos with the same cosmological parameters. 

much larger than the hubble horizon scale, where the effective mass is defined by m^^ — • 
To remedy this situation one should consider an appropriate quintessential potential which 
has a mass comparable the horizon scale at present, and the models considered in this paper 



are the case ISj. Interestingly, some authors have found that one can construct viable 
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MaVaNs models by choosing certain couplings and/or quintessential potentials [48|, 
Some of these models even realises mefj ^ H . In Fig. (fTTl) . masses of the scalar field relative 
to the horizon scale rries/H are plotted. We find that rrics < H for almost all period and 
the models are stable. We also dipict in Fig. (fTT|) the sound speed of neutrinos defined by 
Cg = 5Pp/5pi, with a wavenumber k = 2.3 x 10"'^ Mpc~^. 

In summary, we investigate dynamics of dark energy in mass- varying neutrinos. We show 
and discuss many interesting aspects of the interacting dark-energy with neutrinos scenario: 
(1) To explain the present cosmological observation data, we don't need to tune the coupling 
parameters between neutrinos and quintessence field, (2) Even with a inverse power law 
potential or exponential type potential which seem to be ruled out from the observation of 
uj value, we can receive that the apparent value of the equation of states can pushed down 
lesser than -1, (3) As a consequence of global fit, the cosmological neutrino mass bound 
beyond ACDM model was first obtained with the value '^rrii, < 0.87e\^(95%CL). More 
detail discussions and theoretical predictions on the equation of state and on the absolute 
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FIG. 11: (Left panel): Typical evolution of the effective mass of the quintessence scalar field relative 
to theHubble scale, for all models considered in this paper. (Right panel): Typical evolution of the 
sound speed of neutrinos = 8Py/5pi, with the wavenumber k = 2.3 x 10^^ Mpc~^, for models as 
indicated. The values stay positive stating from 1/3 (relativistic) and neutrinos are stable against 
the density fluctuation. 



mass bound of neutrinos from beta decays and cosmological constrains will appear in the 
separated paper |5l[ |. 



APPENDIX A: CONSISTENCY CHECK 



The form of /t can be also obtained by demanding conservations of energy and momentum, 

(0) 

i.e., demanding that STj^+W ^ STj^ = 0. Let us begin by considering the divergence of 
the perturbed stress-energy tensor for the scalar field, 

( ■■ ■ (PV 



(Al) 

where in the last line we used eqs.(l3]) and (|T9l) . The divergence of the perturbed stress-energy 
tensor for the neutrinos is given by. 



6T^= ~6p - (p + P)d^Vi - m{5p + 5P) - ]-h{p + P) 



for the time component and 



5Ti:= (p + P)v, + (p + P)v, + 4H(p + P)v, + d,P + d,^ 



(A2) 



(A3) 
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for the spatial component. Let us check the energy flux conservation for example, starting 
with the energy flux in neutrinos (in fc-space): 



(A4) 



where 6y = ik'^Vy j. Differentiate with respect to r, we obtain, 

{pu + Pu)eu + {pu + Pu)eu = ^Trka-^ I q^dqqfo^ 



-An{p, + P,)d, (A5) 
Let us consider the flrst term in the right hand side of the above equation. This gives 
Anka"^ j q^dqqfo^i 
= Anka"'^ / q^dqqf^ 



3 e 



kHP, - k\p, + P,)<J, + ^Tik'^a-^ / q^dq^^ 



+ Aiika ^ / q^dqqfoK, 



q^ de 
e 



where a is deflned as (p + P)cr = —{kikj — and expressed by the distribution function 



as 



{Pu + Pu)(^p 



in 



q^dq^h{q)^>2 



(A6) 



Comparing eq. flASP with eq. flASp . we flnd that the divergence of the perturbed stress-energy 
tensor in spatial part for the neutrinos leads to 



d'V^ 6X1"= lAnk 



2 -A 



q^ da 

(fdq——5(j)fQ + A'nka~'^ I q^dqqfoK 



(A7) 



3 J ' 

On the other hand, the divergence of the perturbed stress-energy tensor in spatial part for 
scalar fleld is, from eq. OAip . 



{Pu - ^Pu) 



-ATTk'^5(f)a ^ I q'^dq^fo . 



(A8) 



These two equations imply that k, shold take the form as eq. fl32]) . 

Next let us check the energy conservation. Density perturbation in neutrino is, (see 
eq.dSHD) 



Spu = a 



(271)' 



re/o(g)*o + a 



(27r)3 50 



5<Pfo 



(A9) 
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By differenciate with respect to r, we obtain 



-4 / (1^(1 -r .T, , -4 f d^Q 



d?q d f de 



(27r)3 9r \d<p 



where 



fnj 2 2 I 2 2^^^"^^ 'l\ I 

e = [rCa m + a m — 0)/e 



^ ( - 1 = -^-^7u + + (All) 



dr \d(j) J e 
Inserting eq. pHj) for \l/o in the above equation, we obtain 



Comparing with eq.( lA2p . we find 



which is found to be equal to -V^ 6Ti;= -6 f ^ j - ^50. 

APPENDIX B: BOLTZMAN EQUATIONS IN INTERACTING DARK 
ENERGY-NEUTRINOS SCENARIO 



From the lagrangian C = m{(f))^y—g^uXfJ,X|^, the Euler-Lagrange equation is given by 

d f dC\ dC 



dX \ dx^^ J dx^^ 



Bl) 



where 



dC 


= P, = - 


dxf^ 


dC 


dm 1 


dx^^ 


~ dxi^y 



X- 



\/-gaf^X"'X^ ' 



(B2) 



-g^px-xP - m(x^) 7"'^ = (B3) 
2^J-gapX°'xP 
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Therefore eq llBip becomes 

^-ga^x-xPdX \ ^-g^px-xP j 2 g^^x-x^ dx'^ 

By using the relation ds = — ^ —gapx^^x^dX^ we obtain 

pf^ = -m(x^) , = mix^") ^ (B5) 

and eq. (lB4p becomes 

d / , ,,. dx^\ mix^) dx"" dx^ dm , 

With simple calculation, finally we obtain the relations: 



ds ds " 

dr 



(J I 

P^-^ + Tl^P''P^ = -mg^f'm,^. (B9) 



For /i = component, eq. (]B9P can be expressed as 



l^{Pr + T% = -mg'^m,.. (BIO) 



Since P° = g^'^Po = a2e, each terms of the eq. (lB10p are given by: 



First term = -2a~'^Hq^ + a~'^q'^ - a'^Hm'^ + a-'^m^ (Bll) 

dr dr 

Second term = 2a-^Hq^ + a-^Hm^ + a-^hi^^q^q^ (B12) 

Third term = a~'^m—— (B13) 

or 

Since the first term includes the total derivative w.r.t. comoving time, we obtain finally the 
eq.(H2|) in Section III-C: 

dq 1 • o-mdm dx^ ^, 

_ = „V-a^-^-. (B14) 
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